Nonlinear PDEs for gap probabilities in random 
matrices and KP theory. 

M. AdleiQ M. CafassdU P. van Moerbeka!] 
Abstract 

Airy and Pearcey-like kernels and generalizations arising in random matrix 
theory are expressed as double integrals of ratios of exponentials, possibly mul- 
tiplied with a rational function. In this work it is shown that such kernels are 
intimately related to wave functions for polynomial (Gel'fand-Dickey) reductions 
or rational reductions of the KP-hierarchy; their Fredholm determinant also satis- 
fies linear PDEs (Virasoro constraints), yielding, in a systematic way, non-linear 
PDEs for the Fredholm determinant of such kernels. Examples include Fredholm 
determinants giving the gap probability of some infinite-dimensional diffusions, 
like the Airy process, with or without outliers, and the Pearcey process, with or 
without inliers. 
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1 Introduction 

The purpose of this paper is to show that, given an interval £cS (or a disjoint 
union of them), the Fredholm determinants of kernels of the type: 



(2ni) 2 J T + J r ~ " J e-^W^ V^~W it-?/ 

where are appropriate contours in the complex plane and V p (u) is a poly- 
nomial in u of degree p + 1, properly parametrized by ti,t 2 , ■ ■ • , satisfies 
some non-linear PDEs, which are generated by KP integrable hierarchies. The 
Fredholm determinants of such kernels describe in some cases the gap probability 
of certain infinite-dimensional diffusions, like the Airy process, with or without 
outliers, and the Pearcey one, with or without inliers. This area has been studied 
by many authors, see for instance E51 12HI ESI El E51 SOI LT1 El HI LI2] in the 
case of the Airy process, or 113 113 QSl 021 EH E3 El El Oil 02] in the case of 
the Pearcey process. 

For example, consider A Brownian motions starting and ending at x = 0, 
going from time t = to t = 1, and conditioned to not intersect. We are interested 
in the case when N is very large. Then the average mean density, at time t e 
[0, 1], has its support on the interval I t := — y/2Nt(l — t), y/2Nt(l — t) , hence 

forming a set bounded by an ellipse in the (x, t)-plane. If we look at any point 
on the ellipse with an 'Airy microscope", i.e. with spatial scaling A -1 / 6 and time 
scaling A -1 / 3 (compatible with Brownian motion), then in the large A limit we 
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find [26, ',<£[] the stationary Airy process A(w), whose gap probability is given by 
the Fredholm determinant 



P (A(w) n E = 0) = det (1 - XE Kf m ) 
(here and below xe represents the indicator function of E) with 

tr r - — +X'u / \ n 1 

1 / . / . e 3 ^ I u — w\ 1 



(2m) 



du / cif 



e s 



v — w 



u — V 



1.2) 



;i.3) 



a special case (p = 2) of (1.1) and are contours passing through the origin 



and extending to infinity with an angle of tt/3, see also Corollary 2.9 below. 
It is a consequence of (1.2) and Theorem |1.1| of this paper, formula (1.21), that 



Q = Q(E) := log¥(A(w) n E = 0) 

satisfies the PDE 

d 4 Q + 6(<9 2 Q) 2 + (2 - Ae)dQ = (1.4) 

where, given a disjoint union of intervals E := Ui=i[ a 2i-i? a 2i] C K, we define the 
two differential operators 



<9a,- 



;i.5) 



Now consider, in our non-intersecting Brownian motion model described above, 
the following modification. Instead of requiring all the particles to end, at t = 1, 
at x = 0, we allow n outliers to end at x = e~ ro a/ N/2 and we look with the Airy 
microscope described above about the point on the ellipse given by 



(a, t„) = \/2Nt Q {l - t c 



2(tq+w) 

1 + e jv-i/3 



This yields, in the l arge iV limit, a new process [3], namely the n-Airy process 
A (n \w). Instead of (1.2) we now find, for the gap probability, 



P (A n (w) n E = 0) = det (1 - X eK^ w) ) 
and using Theorem [TTj JT25| ) with p = 2, we find that 

Q = Q(E) := \og¥(A n (w) f] E = 0) 

satisfies the PDE 

d 4 Q + 6(<9 2 Q) 2 + (2 - 4(e - io^))5Q + 3<9 2 Q 



1.6) 



0. 



;i-7) 
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Both the equations (1.4) and (1.7) are new (here and in the paper we denote, 
given an arbitrary variable w, d w := ^). In [5] another much more complicated 
but nevertheless useful PDE for O was found. 



Another beautiful example, again in the framework of N non-intersecting 
Brownian motions, involves all particles starting at t = at x — 0, with half of 
them ending at -y/N/2, the other half at y/N/2. Now for a time t e [0, 1/2] the 
average mean density is supported on one interval that, for t > 1/2, separates into 
two. In the (x, t)-plane the edge of the support is heart-shaped, with a cusp at 
(x, t) = (0, 1/2) and the heart is anchored at the points (0, 0), (±yjN/2, 1). If we 
now look about the cusp at (0, 1/2) with a "Pearcey microscope", i.e. with spatial 
scaling iV -1 / 4 and time scaling iV -1 / 2 (compatible with Brownian motions), then 
in the large N limit we find [37] the Pearcey process V{r) with gap probability 
given again by a Fredholm determinant (here E must be compact) 



p (V(r) n E = 0) = det (1 - XE KT m ) 



'U 



with 



K^ w) {\>:):- — 



4 2 

C 4+ T 2 / U + W 



du I dv 



u — v 



1.9) 



another special case (p = 3) of ( 1.1[ ) (see Corollary 2.9 for a description of the 
contours). In this case we get for 

Q = Q(E) := logP(P(r) n E = 0) 
the following set of equations, obtained from (1.8) and Theorem |1.1| upon spe- 



cializing equations (jl.21 )-( 1.24) for p = 3 and upon setting r 

<9 4 Q + 6(<9 2 Q) 2 + 12<9 2 Q - 4r«9 2 Q = 



2U 



2d T d 3 Q + \2{d T dQ){d 2 Q - r/3) + (1 - 3e + 2rd T ) 

8d 3 T Q + 4 {d T dQ, d 2 Q} 9 + (e - 2rd T - 2)d 2 Q = 

d A U + \2dlU + 6d 2 U 2 = 0, U:=d 2 Q-^. 



;i.io) 

;i.ii) 
;i.i2) 

;i.i3) 



The last one is just an other form (Boussinesq form) of (1.10). Equations (1.10) 



and (1.11 ) are new, while (1.12) was first derived in [7j. In the above and through- 
out this paper, we denote {/, g} g := f(dg) — (df)g. 



x 



In the previous Pearcey example, if we allow n inliers to end, for t = 1, at 
= — w (y) 1 ^ 4 , that would lead, in the large iV limit and near the cusp, to 
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a different process, namely (see [3]) the n-Pearcey process pw(r). The gap 
probability of this process (again E is compact) is given by 

P (?W(t) n J5 = 0) = det (1 - ^ )XE ) , (1-14) 
leading by Theorem exactly as for Pearcey, to the equations ( 1.26 )-( 1.29) for 

:= logP(P {n) (r) HE = 0), 



upon setting r = 2t%. All of these four PDEs are new, although the third one 
was conjectured in [6]. 

It is also worthwhile to mention that the PDEs that have been derived can be 



used to obtain asymptotic results. For example, as a consequence of (1.12), one 



can show in the style of pQ the following approximation of the Pearcey process 
by the Airy process, by looking far out along the Pearcey cusp to yield, for large 



r: 



p i t(t) - mi n { _ E) 

(3r)s 



P(^(0)n {—E) = 0) (l + 0(r 



(1.15) 

where as usual E is a finite union of compact intervals (see also [T2] for a similar re- 
sult obtained using the non-linear steepest descent method for Riemann-Hilbert 
problems) . 

In order to state the main theorem of this paper with precision, we need the 
following preliminaries. 

Given the root of unity uj = e^ 1 , consider any subset of rays taken from 
configurations 

-p + 1 



J 



< 



}, r p c{c(u^), j< [f]}, (i. 



16) 



consisting of oriented counter-clockwise contours C(u 
integer n > 0, consider t 

41w( A ' A, ) = ^ I du I dv 



an integer n > 0, consider the kernel in (j 1 . 1 p 
where V p (u) is a polynomial 



= + Given 

e -V p (u)+X'u f u _ w \ n 1 



Vp(v) + \V \ y 



W U — V 



(1.17) 



p+l i+l 

V p (u) = ^—7 + Y J d l ^—, 6i :=^(t 1; ...,Vi: 
p+l ^ I + 1 



P+i 



p-3 



(*1 



'1.1? 
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with polynomial coefficients 9i in t implicitly given in terms of ti,...,t p _i, by 
solving the equation w = Vp(u) for u in terms of a series in large w, as in (i) and 
identifying with another series, as in (ii); thus 

u = W p 6> p _ 2 w p Opsw p — 9 2 + 9 p -4 w p 

p p p \ 2p y J 

+ ■■■ + 0(w~ 1 ~p) (1.19) 

= wp + - J2(p - j)tp-jw~p + o(«r 1 - 1/p ). 

V 1 

The Hirota bilinear equations for the KP-flow imply two strings of partial 



differential equations, characterizing the KP r- function, namely (see Lemma 4.1 
Y<? : (pe+i(d t ) - -d x d^To t = 0, Yi^-i : [d x d t - ]-d 2 d^ x - d x pi{d t 



t o r = 0, 



We further remind the reader of the definition (|1.5|) and the previously stated 
Theorem 1.1 Each of the Hirota equations 



convention that di := -Sr and d w = Now we can finally state our main result. 



Y 3 , . . . , Y p+1 , Y M , . . . , Y ljP and 2{p + 2)Y p+2 + (p + 1)Y 1)P+1 , 



gives rise to a non-linear PDE for the Fredholm determinant of the kernel (1.1) 

Q = Qpfo, Vil E ) ■= logdet(/ - ^....^^l^o- (1-20) 

These PDEs only involve the differentials d and e with regard to the boundary 
points of E and the t-partials 82,83, ...,d p - X . Some noteworthy examples are: 



lr £Yie Hirota symbol of two functions / and g, associated with any polynomial pit), is given 

by 

while the elementary Schur polynomials are defined by the relation 

00 
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Case 1 (n = 0) 



• For all p > 2, Q satisfies 
the Y 3 -equation: 

d A Q+6(d 2 Q) 2 + 5 2 , p (2- 

+ {l-S 2tP ) (zdlQ-A 



3- — -t p -id + (1 - 5 3 , P )<9 3 
V 



(1.21) 



• For all p > 3, Q satisfies the equation (1.21) and 
the Y 4- equation: 



P 



{p - 2)t p . 2 d + (1 - 5 A , p )d A 



(1.22) 



+(1 - <f 3sP ) [2d 2 dzQ - 3 

the d 2 Y 3 — dY \- equation 

d 3 2 Q + 2 {d 2 dQ, d 2 Q} 9 - 2(1 - 5 3>p )d 2 d 3 dQ + 5 3jP (e - 2t 2 d 2 - 

2 ( v \ o (1-23) 

- -(1 - <W) ((p - l)Vi9 2 - 2(p - 2)V 2 9 - |(1 - 5 4 , P )<9J <9 2 Q = 



and the d 2 Y 3 - equation (the Boussinesq form of the Y 3 -equation) 



d*U + 3d 2 U + 6d 2 U 2 - 4(1 - 5 3 , p )d 3 dU = 0, where U := <9 2 



p — 1 
P 



Vi- (1-24) 



Case 2 (n > 0) : Q satisfies 



For p = 2: 



For p = 3: 



<9 4 Q + 6(d 2 Q) 2 + (2 - 4(e - wd w ))dQ + 3d 2 w Q = 



<9 4 Q + 



8t 2 d 2 Q + 3<9 2 2 Q + 4a, 



0. 



(1.25) 



(1.26) 



(d 2 d 2 - 2t 2 d 2 - 3(e - wd w ) + 1)<9Q + 6(<9 2 Q)(<9 2 <9Q) - 2Q 2 U ,Q = 0, (1.27) 



(e - wd m - 2t 2 d 2 



+ d 3 2 Q + 2{«9 2 dQ, d 2 Q} 9 + 2d 2 d v 



and <9 4 £7 + 3d 2 U + 6d 2 U 2 + Ad w dU = 0, 



U = d 2 Q - §t 2 . 



0, (1.28) 
(1.29) 



Notice that for n = 0, Q a/so satisfies these equations, but without the terms 
containing d w (as it does not depend on w). 
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Theorem 1.1 is new, and it has as consequences the PDEs (1.4), (1.7), (1.10)- 
(1.13), and ( |1.26 )-( 1.29 ), for the Airy and Pearcey cases, with or without in- 
here and outliers, as previously discussed. However, at heart, the theorem is a 



method for yielding PDEs for Fredholm determinants (1.20), with kernels of the 



form (1.17), some of which were derived explicitly for the benefit of the reader; 



although clearly our motivation is random matrix theory. Indeed it would be in- 



teresting to relate all the kernels in (1.17), also for p > 4, to the theory of random 



matrices (see open questions below). Other methods for computing PDEs for gap 
probabilities can be found in [3 El [3 El El EH [36l HH E21 H21 [HI [33l [32] . 

The main tool of this paper is the natural integrable deformations of the 



kernel (1.1), which enables us to represent it essentially as an x-integral of the 



product of a KP-hierarchy wave function and its adjoint, and hence in terms of 
KP tau functions, as elaborated in Proposition |3.3[ Consequently, we can use 



the the tools of Sato's (formal) KP theory of Grassmannians, which are provided 
for the reader in Section 2, along with applications of these tools, which we shall 
need in the next section. In Section 3, we show, using a proposition originally 
appearing in how to use the fact that the "integrable deformation" of the 



Fredholm determinant (1.20) has an expression as a ratio of two tau functions 
as in (3.8), (in essence a continuous soliton formula), to derive the so-called 
"Virasoro relations" for the Fredholm determinant itself. This leads in Section 

or more precisely to a method to 



4 to a derivation of the PDEs of Theorem 1.1 



derive PDEs. In the appendix we apply the method to derive the first example in 



Theorem 1.1 , leaving the others to the interested reader. There is a large number 
of formulas in this paper, first going from Sato's theory to its applications, and 
then in assembling the formulas necessary to derive any PDE inherent in Theorem 



1.1 Removing some of the formulas would make it impossible to check the results 
of the paper and moreover to reproduce these PDEs. 
We end this introduction with some open questions: 

1. In this paper we derived PDEs just for scalar Fredholm determinants re- 
lated to gap probabilities, not for the matrix ones describing joint gap 
probabilities for multi-time processes. It is an open problem to develop the 
integrable theory for the matrix Fredholm case and hence derive PDEs for 
multi-time joint gap probabilities (see also |llj . where it has been proven 
that the multi-time Airy and Pearcey kernels are integrable kernels in the 
sense of Its-Izergin-Korepin-Slavnov [21]). An example of such PDEs, to be 
derived for the Pearcey process with inliers, can be found in [6], Conjecture 
1.5. 



2. Since the PDEs derived ultimately come from the KP hierarchy, one might 
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expect that they have a "Painleve property" in a suitably general sense; 
for example, the solutions should have only poles along a movable divisor 
in the space of parameters. A first step in this study has been done in 
[12]. Indeed, in that paper, the equations fll.10 ) and (1.11) have been 



re-derived (together with an additional one) as compatibility conditions 
of an isomonodromic system of ODEs. It is a beautiful open question to 
develop a theory establishing a general connection between the equations 
we derived here, the theory of isomonodromic deformations and related 
Painleve properties. 

The question of appropriate initial conditions for such PDEs is also an 
open question and might involve asymptotic solutions rather then analytic 
solutions of such PDEs. This question is not unrelated to the second open 
question. 

According to a well-established physics literature (see [10] and references 
therein) a large class of kernels (the so-called (p, g)-kernels) should appear 
considering critical phenomena in two-matrix models. These kernels should 
be related to Gel'fand-Dickey hierarchies supplemented with appropriate 
Virasoro constraints, even if no general expression is given yet (see also [19] 
and references therein for the (2, g)-case). It would be extremely interesting 



to prove in full generality that the kernels (1.1) we studied here, given an 
appropriate choice of the contours, correspond indeed to the (p, l)-kernels, 
as it is the case for p = 2, 3. We also mention that, using the method 
developed in this paper, in [2] we have been able to derive PDEs for some 



kernels generalizing (1.1) and including the ones studied in [19] . 



2 KP Theory revisited 

In this section we provide the basic tools, which will be needed, (see for example 
[20] |2T] HI |3l [29]) from KP theory; in particular r-function's formalism and wave 
functions, Grassmannian theory and additional symmetries, and also p-reduced 
(Gelfand-Dickey) and fc-vector p-reduced KP hierarchies. These latter intro- 
duced by Krichever [2B] and more specifically we will use some results due to 
Helminck and van de Leur (23] [38]. We shall just give the basic facts necessary 
for the subsequent chapters. 

Let q = q(x, y, t) be a function of three variables; the KP equation describing 
shallow water waves in M 2 , namely 

q xxxx + I2ql + I2qq xx + 3q yy - Aq xt = 0, 
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is actually part of a hierarchy of commuting PDE flows for 

q := dl\ogT(t 1 ,t 2 ,t 3 ,t 4 , . . .) 

where we should identify^] 

ti =x,t2 = y,t 3 = t 

and the other "times" corresponds to higher order flows of the hierarchy. The 
whole set of equations satisfied by the hierarchy can be written, using the function 
r(t) := r(ti, t 2 , ta, t±, . . .), as the famous bilinear identity 



r(t - [z- 1 ])^' + [z-^fPrw = 0. 



(2.1) 



In the expression above we denoted [z] := (z, z 2 /2, z 3 /3, . . .) (and similarly for 
z^ 1 ); the equations of the hierarchy are obtained expanding the integrand as a 
formal Laurent series about z~ l = and then taking the (formal) residue about 
z = oo. 

The hierarchy can also be written as a commuting set of Lax equations for a 
(x, t)-evolving pseudo-differential operator 



L + = L + (x, t):=D + ^2 afa, t)D~ 

3=1 



(2.2) 



where D := d x (the reason for the notation L + instead of the classic L will be 
clear later) ; namely 



d_ 
dti 



i e Z , 



(2.3) 



where, for an arbitrary (formal) pseudo-differential operator H := Xliez hiD 1 , we 
denote 

{H) + :=J2hiD\ H := H — H + 

i>0 



d 

(note also that from the equation (|2.3|) we deduce = D). 



2 actually in the sequel we prefer to keep t\ and x as two distinct variables; the identification 
is justified since r always depends just on the sum x + 1%. 
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In order to derive the Lax formulation from the bilinear identity (2.1) it is 



necessary to define a (suitably normalized) eigenfunction of L + , itself defined 
by a pseudo-differential operator W introduced below. Indeed setting 

, n ( d 1 d 1 d \ , AS 
t := t + xei = (x + tl , t 2 , t 3 , . .), ft : = - .j (2.4) 

and denoting as before with Pi(t) the classical Schur polynomials such that 



^=^^(t), 

i=0 



we define the wave operator W = W(t) 

w l= Tjt-jD- 1 }) ^^ _ ( ^ Pj (-d t )r(t) 



r(t) 



E- 



Tt 



(2.5) 



(2.6) 



Denoting by if* the formal adjoint of a given pseudo-differential operator 
(the formal adjoint acts through the formula (a(x)Di)* := (— DYa(x), j G Z 
and linearity), we can define the wave function and Lax operators L + ,M + 
together with their adjoints (denoted with the minus sign) 



L+ ■= WDW- 



M+ := \VxW~ 1 . 



(W 



(2.7) 



— i \* / 



-Z))W*, M" := (W^YxW*, := (t^" 1 ^"* 2 



One should think of L 1 * 1 and M 1 * 1 as the "dressing" of ±D and x while the wave 
functions ^ are the "dressing" of e ±:C2 ; this leading to the following relations, 
which result from "undressing" the operators and functions, to wit: 



±iT/± 



d 

z y± M ±ty± = ± q± \ L ± M ±] = ±1 . 
oz 



(2.8) 



Also the definition of the wave operator W (eq. (2.6)) together with (2.7) gives 



back the celebrated Sato's formula expressing the wave functions in terms of the 
tau function, namely 

p ±xz±y:Zi Uz* r(t =f K 1 ]) 

Til) ■ 



(2.9) 



It is a well known result that Sato's equation (2.9), together with the bilinear 
identity (2.1), gives the following deformation equations for the wave functions^] 

= ±(L ±l ) + v|/±, ieZ+ , (2.10) 



3 Here and below L 1 means (L ) l and not L to the power —i (on the other hand we did 
not introduce any operator L). 
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and these equations, finally, give as compatibility conditions the Lax equations: 

Al±=[±(l±*) + ,L ± ] *GZ + . (2.11) 

As a side remark we want to point out that, till now, the pseudo-differential 
operators M did not play any role. As a matter of fact they are not useful to 
define the hierarchy, but rather to define some additional symmetries which will 
be discussed later. 



Sato observed that the KP flows linearize on an infinite-dimensional Grass- 
mannian Gr whose elements are linear spaces W of formal Laurent series, i.e. 
typically we have^] 

W = span I ip n (z) = z n + a ni z\a ni ec\ . (2.12) 

L -oo<i<n— 1 J n=0 

More precisely, associating to any KP solution, through its wave functions, two 
subspaces W ± through the formula 

W ± := spaB i > {D i * ± (a;,0;z)}, (2.13) 

we get that the linear spaces W are t-deformed by the KP flows via 

W ± (t)=e^^ t ^W ± . 

In fact the subspaces W ± defined above are an alternative way of characterizing 
the KP data encoded in the tau function (or also in the wave function). In 
particular invariance properties of the KP solutions are sometimes more easily 
characterized in terms of W ± . To this end consider the mapping from z-operators 
A to x-operators Va given by 



-oo<i<oo j>0 v / i,j 



— oo<i<oo j>0 ' i,j 

(2.14) 



4 As the expert readers will have noticed, we are just describing the so-called "big cell" of 
the Sato's Grassmannian, which is enough for our purposes. 
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where the equalities are proved using (2.8). The following identical statements 
follow from (2.14) and relate the invariance properties of W ± with regard to 
differential operators in z to properties of the associated differential operator in 



A+W + C W 



V 



A+ 



(V 



A~W c W 



(2.15) 



z 2 are 



In particular the points in the Grassmmannian satisfying (2.15) for A A 
the points associated to the KdV hierarchy, considered as a reduction of the KP 
one (more details will be given in the following section). It should be noted that 
Sato's Grassmannian, differently from the Segal- Wilson one |34j . is constructed 
with formal series, not analytical or asymptotics ones. Hence all the relevant 
computations are done with formal series, and in particular the basis in (2.12) 
are well defined as formal series while the tau function r(t) is a formal infinite 
series composed of monomials in the ti, computable from the formal basis; in 
the sequel we will take full advantage of that. Nevertheless, for some particular 
subspaces W and some particular values of time parameters, formal computations 
acquire a more concrete meaning as analytical and/or asymptotics series, (on the 
other hand it is worth recalling that the Segal- Wilson Grassmannian is strictly 
contained in the Sato one). 



2.1 p— reduced KP— hierarchies 

In this section we recall the definition (in this specific context of Sato's Grass- 
mannian theory) of Gel'fand-Dickey hierarchies which we shall refer to as the 
p-reduced KP hierarchies. First we need the following well known lemma: 

Lemma 2.1 Let p be a non-negative integer. Given a point W + in the Sato's 
Grassmanian and let L + and r be the Lax operator and the tau function of the cor- 
responding solution of the KP hierarchy. The following conditions are equivalent 
and are conserved along the KP flows. 

1. z p W + C W+ 

2. L +p = (L +p ) + 

3. t does not depend on t np for any n > 1, modulo a removable factor of the 
form e^^ 1 c " tnp . 



Proof: The equivalence between 1. and 2. is given by (2.15 ) while the one between 
2. and 3. can be proven observing that, if L +p = (L +1 



then we have 



d 



Ot 



Vn > 1 



np 
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and so 3. is equivalent to 2. using Sato's formula (2.9). The fact that the three 
conditions are conserved along KP flows is proven by observing that the first one 
is obviously conserved along the KP flows, since W + (t) = e"^ 1 '' 2 '^ in Sato's 
Grassmannian. ■ 

Definition 2.2 The p-reduced KP -hierarchy is the KP-hierarchy supplemented 



with one of the (equivalent) conditions 1. 2. or 3. of Lemma 2.1 



In the case of p-reduced KP hierarchies the relevant Lax operator, rather then 
being L + , is its p th power L +p , since the latter is a differential operator. Indeed, 
upon using the Lax formulation of the KP hierarchy and setting 

:= L ±p , 

one finds equations 

£ ±^± = z Pty± } = (±C ±l/p ) + V, i G Z+ 

Obi 

and their compatibility conditions give the Lax formulation of the p-reduced KP 
hierarchy: 

f)r± 

— = [{±£±"*) + ,£], *GZ + . (2.16) 

In the sequel we are interested in some particular solutions of p-reduced KP 
hierarchies satisfying an additional invariance property. Let's start defining a 
^-operator 

it is easy to check that Ap satisfies the following condition: 

iA±(z),z p ] = ±l. 

Theorem 2.3 The invariance conditions 

z p W + C W + , A+W + C W + , (2.18) 

determine uniquely a plane G Gr, which moreover uniquely determines W~ G 
Gr by the relations 

z p W- cr, A p W~ C W~. (2.19) 

The Wp are linearly generated by the eig en functions cp^ of the operators Ap, 
namely: 

W± = span,> {(^J)Vj} (2.20) 
13 



with 



OO -|- 



(2.21) 



the latter uniquely determines (Pp(z). 

The corresponding p-reduced KP wave functions ^p(x,0;z) are then uniquely 



specified by 



*±(0,0; = 



d 



(2.22) 



A±{z)*±{x, 0; z) = ±^ 0; z). 



Parts of this theorem appeared in [27] and jl], but perhaps with sketchy proofs 
for the case p > 2. We will see in the sequel that these particular solutions of 
KP, determined either by Wp , y?p or ^ (a;, 0; z), can be made fairly explicit (see 
Theorems 2.6|2.8 and Lemma 3.7). 



To prove the theorem above it clearly suffices to prove the following lemmas. 
Lemma 2.4 The condition 

z p W + C W + , A+W + C W + , (2.23) 

uniquely determines W + together with its adjoint W~ , and these are specified 
by (2. 2(fy, w ith ^fp(z) determined uniquely by the differential equations and the 
series \2.2~ty. 



Lemma 2.5 The wave functions ^p(x, 0;z) going with Wjf of ggp satisfy the 
equations (2.22) and are uniquely determined by them. 



Proof of Lemma \2.4 ' Note that since 



At 



z + O 



W + e Gr, 



then for some <p+ (z) = 1 + we have that 

i=i z 



W 



span { 04J") V P 



O 



but since z p W + C W + , z v Lp J p r (z) = z p (l + O (^)), we must have, for some con- 
stants Co, ci, c p _2, that 



z*v; 



( W + c^+r 1 + Cp _ 2 (^+)^ 2 + ... + c ) <p; 
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and so we just must show all the Cj = 0. Sincejj A£ = z + by induction 

(At) p ~ j = z p - j + e-^-i, and so 



= ((^ + Cp _ 1 (^-i + Cp _ 2 (^+)^ + ... + Co _ z P)^ 



(£_! + c p _!(2; p 1 + £_ 2 ) + c p „ 2 {z p * + e_ 3 ) + ... + c x [z + e_ p ) + c ) 



y p-2 



So we first conclude c 



p— i 



0, since the left hand side has no matching power 



of z p 1 , then c p _2 = 0, and so inductively, cq = 0. Finally (2.15), (2.23) implies 
z p W~ C W~, A p W~ C W~ and we can repeat the above arguments for W - , 



concluding the proof of Lemma 2.4, except that it remains to be shown that (2.21 ) 



determines <f^ (z) uniquely. To do so, from the above (Ap) p = z p + e_i, and more 
precisely 

{{a;y -z p ) = j- + 5_ 2 + <l 3 + ... + <*_ p2 , 

where 5k is an operator such that deg(o~fc(z 1 )) — i + k, for — oo < i < oo. We find 
from ( pE2~T| 

= {(At) v -z p )^{z)- 



d_ 

dz 



+ 5-2 + 5-3 + • • • 



1 + 



+ 



(-S + mu) + (-* + ^(*) + mi: 



+ 



3a, 



6. 



+ 5- 



+ 5_ 



+ ...=: r_ 2 + r_ 3 + r_ 4 + . 



zr \ z / \ z 

with degree Ti = i. Since all the terms Tj individually must be 0, we recursively 
determine a1~ from r.,_3 = 0, j > 1, and similarly for (z). ■ 



Proof of Lemma 2.5: From G Gr, and comparing (2.20) and (2.23), 

namely sp&n i>0 {Lr$T~(x, 0; z)} = span j>0 {(. Ap ) V p }, and observing both (i) 
and #+(0, 0; z) are (l + O Q)) (from (|S),( |2j) an d r(x 0,0 . .) ^ 0) conclude 
^+(0,0; z) p = <p+(z). Moreover, from ( 2.14[ ), ( |2.23~] ) and ( |2.15[ ) conclude, in par- 
ticular, that 



a;(z)v;(x,o;z) 



(V^+(x,t;z)) 



t=o 



(P++(x,0;z)) + ^+(x,0;z). 



'Here £&(•) stands for a differential operator having the property that £fc(z l ) 



-00 < i < 00, the af' being constants 
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Now denote with , M+ the two Lax operators corresponding to this particular 
solution of the p-reduced KP-hierarchy; from (2.14) and (2.17) we find 

1 



and thus 



it=o 







p 



-)Wr p 



lt=0 



p '+lt=o 



d_ 

dx 

(2.24) 



A;(z)*;(x,0;z) = — *;(x,0;z), #+(0,0; z) = ip+{z), 

the latter being a nonsingular first order PDE with given initial condition, which 
completely determines *&p(x, 0; z), as claimed. Meanwhile from (2.14), (2.15), 
fl2.23p and ( |2.17[ ), conclude that 

A p (z)^(x,0;z) = (p-_*;(x,t;z)) 



t=o 



and then as before # (0, 0; z) = y? p (2;), with # (a;, 0; z) uniquely determined as 



before, concluding the proof of Lemma 2.5 



2.2 Airy-like integrals and asymptotics 



In this section we show how to explicitly construct (in Theorem 2.6) the wave 

Their asymptotics is given in 



2.3 



functions j3/^(x, 0; z) described in Theorem 
Theorem 2. 8| together with some examples in Corollary 2.9 
As an application of Theorem 2J3 consider the functions 

yP+ 1 



±p 
2^ 



±uy d 



y, 



(2.25) 



with Tp being a union of lines in C through the origin, picked so that (2.25 ) makes 
sense, and such that there exists regions T>+ which are unions of sectors about 
the origin in which $+ have the following asymptotic behavior (see Theorem 2.8) 
in T>\ 



±r z p ) 



2 QP+1 



afz~ 



for z G Z>±. 



(2.26) 



Theorem 2.6 T/ie wave functions ^>p(x,0;z) of Theorem 2.3 are given by the 
following formula: 

typ(x,0;z) = z^e+p+i 
Moreover they satisfy the spectral equation 

£j(x, 0)*J(x, 0; z) := ((±£>)" - x)*±(x, 0; z 



* P+ V(z + z p ). 



^(x,0;z). 



(2.27) 
(2.28) 
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Proof: Using Theorem 2.3 we simply have to prove equations (2.21) and (2.22). 
Note that by our choice of , c p := \J ±2n/p, the functions $p are solutions of 
the spectral problem: 



=F ± 



d 

du 



± v) $J(u) = - / (TV P ± u)e T i^ ±uy dy 
J c p ir± 



d_ 
dy 



:3 £±^ j rf y = 0, 



(2.29) 



which implies ( |2.28[ ). Note that ( ]2.26[) s hows that ¥±(0,0; z) = 1 + JXi of 
and so to prove the first relation in (2.22 ) we need to prove the following equation: 

{A$y{z)*$(0, 0; z) = z^±(0, 0; *). (2.30) 
To that end, observe, using (2.17) and (2.26) that 
d 



d^> 



= z p pzP- 1 dz py ' 



d 



and so conclude that 
d 



Z 2 



u=zP 



pz p_1 dz 
z - E T i e ± ^ zP+ Uj(^)^(0,0;z), 



1 e ± ^^ +1 (^(z)^(0,0;z)). (2.31) 



Repeating the argument p-times and using (|2.26|) and (|2.29|) conclude that 
-( E r i ) e ± ^ +1 vl/±(0,0;; 



<(«)) 



u=zP 



±- 



du 



u=zp 



((^rw^(o,o^)) 



and this yields (2.30), which yields (2.21) and the first equation in (2.22), while 
the second equation is a consequence of the calculation (2.31) with u y x + z p 
instead of z p . Theorem 2.8 yields (2.26), which was used to derive the series 
expansion for \Pp (0, 0, z) = <fp(z). m 

Prior to proving that the formal series in (2.21) are actually asymptotic ex- 
pansions, we remind the reader of a well known theorem concerning asymptotic 
expansion of solutions of holomorphic differential equations (see Theorem 19.1 of 
[39]). here slightly rephrased for our applications. For C a constant matrix (not 
necessarily diagonalizable) written in Jordan form as C = UDU^ 1 , we denote as 
usual z c := e cl ° sz = Uz D U.~ x and thus z c can be multivalued, depending on D. 
We now have 
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Theorem 2.7 (Wasow) Let A(z) be an n x n matrix function, holomorphic for 
\z\ > \zq\> z a S, where S is an open sector with vertex at the origin. Assume 
that A{z) possesses an asymptotic series "in powers of z~ l ", i.e. A{z) = '^^A i z t . 

i<£ 

Then corresponding to every sufficiently narrow open subsector S , the ODE: 
Y' = A(z)Y possesses a fundamental matrix solution in S of the form 

Y ^ = fe(^i) ^diag( e »^'),...,e".^')), 

with p a positive integer, the Pi polynomials in their arguments, C a constant ma- 
trix, and the sum is an asymptotic series. The solution depends on the arbitrary 
choice of the branch of z Y l p . 



The following theorem proves that the formal expansions of (2.21 ) are actually 
asymptotic expansions. It is worth noting that the expansion are independent of 
the sector they occur in. 

Theorem 2.8 (Asymptotics for general p) Consider the integrals 



$±(tO = A /±|- le^^dy (2.32) 



over contours (see Figure^ T^, with u = ep +1 , 

r+ c {c(^), j < p±I] } , r; c {c(u^), j < [|] } . (2.33) 

Then $^ have expansions (independent from the sector) of the form 

<S>±(zn = z- (E ^ ) e ± ^ zP+1 [1+ 4) ( 2 - 34 ) 



i>[m 



z- 

_ 

2 



for sectors of size ^^pjy centered about the rays 



c(u 2j ) c r+ c(J 2 ^) c r- (2.35) 



The contours are all oriented counter-clockwise with the contour 
C(u e ) = R + ■ u e + M+ • tu e as a point-set. 
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Proof: Considering a specific contour C(u 2 ^) C r+ and the sector of angle size 
p(p^) centered about this ray; the points along a ray within this sector can be 
parametrized by u^e z, with z G M + , and |5| < , ^ . Then, setting y i-> yz, 



z G 1R + in the integrand in (2.32), one obtains 

$+((aWz) p ) = ./-^z / e v{y)zV+1 dy, with F(y) = + e ip6 oj 2jp y. 

V 2vr J Fp p+1 

(2.36) 

Employing the saddle point method, first note 

V'{y ) =0^y = e ie u* ■ {1, 7 , ■■■,1 P ~ 1 }, 7 = eT\ 
Then, using u = ep +1 , 

and since 2 G M + , 

3tf(\/(y )^ +1 ) = -^z p+1 cos ( (p + 1)5 + — ] , < £ < p - 1. 
p+1 V P / 

Except for the prefactor ^^z p+1 , the values of $l(V(yo)z p+1 ) are given by 

cos ((p + 1)5) , cos ^(p + 1)5 + , . . . , cos ^(p + 1)5 + 2?r ^~ - j , (2.37) 

with the maximum given by the first number, namely cos ((p + 1)5); that is for 
I = 0. Indeed, when 5 > 0, then cos((p + 1)5) > 0, since < (p + 1)5 < vr/2. 



From there on, the sequence (2.37) goes down and then goes up again, finally up 
to 

cos f (p + 1)5 + 2n ( p ~ LM = cos ( _ (p + i)^ ) < cos ((p + i)^) ; 



p J \ p 

since 7r > 2ii/p — (p + l)5 > (p + l)5, upon using |5| < p ^ +1 ^ ■ A similar argument 
holds when 5 < 0. Thus only yo = e^u 2 ^ , among the p roots of V'(yo) = 0, will 
count in the saddle point analysis, since this point belongs to the sector centered 
about u 2j . So the contour C{uj 2 ^) can be deformed to pick up the saddle point 
yo = e u; 2 - 7 . The function V(y) has the following form, near yo: 



V(y) = V(y ) ~\{{y- yo)V-V"(y )) 2 + 0(y - y ) 3 , with - V"(y ) 



p-i 
VVo , 
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which leads to the new integration variable u G M: 



y = yo + 



'-V"{y G ))~ 1 / 2 r.-V2 1( ( l -f)/ 2 



p+1 

Z 2 



-u = y + 



p ' y 



p+i 

Z 2 



-u. 



So localizing about y , compute, setting z := e l6 ui 2 ^z 



$+(F) = <5>+({u 2l e ie z) p ) 




r p+i (p- 1 ) / / i 

/ e -|(^v^^(^o)) 2 h+oj^ 

ir p V \z 2 



2 2 



_P e« W «)H-i ^/27T ( 1 + (J_ 

eV J ) 1 2 ; e p+ lV ; ( l + O 



£+1 

Z 2 



z ^'ep+i" l + O 



a E±_L 

Z 2 



(2.38) 



Finally to prove (2.34), we must appeal to Wasow's theorem , with our p being 



identified with the p of the theorem. Indeed by (2.29) we find 



y(z) = ($+(*), (*+)'(*), ($+)^(^)) T , r = A(*)y 



1 



o 



o 



1 





yielding (2.34) by Wasow's Theorem and (2.38). 



To do the case of <& , a ssume C(w 2jf+1 ) C T p and in (2.32) compute instead 



$~((u 2j+1 e ie z) p ), replacing (2.36). The function V(y) now reads 



y 



■p+i 



Then one checks: 



e ipd {u 2j+1 ) p y, with V'(y Q ) = for y = e i8 u 2j+1 {l, 7, ...,7 P " 1 } 



P + 



T (* 0^(1, 7> .., 7^} 
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Figure 1: A possible choice of the integration contours for p = 2, 3,4. 



and analogously the root y = e u 2 ^ +1 will now dominate in the saddle point 

analysis. Since T~ D C(lj 2: * +1 ), we can deform r~ to pickup the saddle point 

yo = e u 2 i +1 . The only difference now with the previous case is tha t, when we 

rj- I I (i-p) 

localize, the Gaussian has a positive sign, hence \ appears in (2.32) and iy 2 

i-p 

(instead of y 2 ) is the direction of steepest descent, concluding the proof. ■ 



Corollary 2.9 (Examples of Theorem 2.8) We can pick $p for 
p = 2, 3, 4 as follows (see Figure [7p: 



7T J T 



$ 3 » 




2tt 



2 
JOO 



e- y3/3+uy dy, 
e~ y4/i+uy dy, 



^y 3 /3-vy 



Jr 



dy (2.39) 




2tt 
2~ 



2 j/ 4 /4-t)i/ 



(2.40) 



e- y5/5+uy dy, ®2(v) = J e y5/5 ~ vy dy (2.41) 



which leads to the following asymptotic behavior of (2.26) 



p-l + P ~v 
Z ~—~ ' 



e^Cl + E^-*) 



for p = 2,3,4 in regions Pf 1 . 



2.3 A;- Vector p-reduced KP-hierarchy 

We shall need the following theory of Helminck and van de Leur [231 EE] gener- 
alizing (in the setting of Grassmannians) some known results for p-reduced KP 
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to the case of the rational reductions of KP introduced by Krichever [2H]- It is 
clear that the conditions below, in the Definition 2.10 , are conserved along the 
KP flows, thus leading to a reduction of the KP hierarchy. 

Definition 2.10 The k-vector p -reduced KP-hierarchy is the KP-hierarchy cor- 
responding to points W in the Grassmannian such that 



a subspace W C W of codimension k exists with z p W C W. 



(2.42) 



It turns out that the solutions of the /c-vector p-reduced KP-hierarchy correspond 
to a Lax operator of a particular form. 

Proposition 2.11 Given a solution of the k-vector p-reduced KP-hierarchy let 
L be its Lax operator. Then the pseudo-differential operators £ := (£ ± ) p , satis- 
fying the Lax equations (2.16), also satisfies 



with the functions qj, Tj flowing according to: 

±(L ±n ) + qf, l<j<k, n>l. 



(2.43) 



dqf 



dtr. 



(2.44) 



The case k = 0, namely W = W, reduces to the standard p-reduced KP hier- 
archy; in the next section we will just use the case k = 1, while in this section 
all the results are stated for the general case. As an application of this theory, 
consider, for given Wi G R, 



2^ 



dy, 



(2.45) 



with Tp picked as in (2.25) such that there exist regions Up, unions of sectors, 
in which $^ have the following asymptotic behavior in Up, with n = Y2i n ii 



y ±n \ 2 J g p+1 2 



1 + E 



oo -L 

af 



for z & U 



± 



(2.46) 



Therefore the functions (generalizing (2.27)) 



*J n (x, 0; z) := z^ n+E ^e T ^ zP+1 ^ n (x + z* 



(2.47) 
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behave, for x = 0, as asymptotic series in z 1 , as 

oo 

*J n (0, 0; z) = 1 + j3 



(2.48) 



Now define the operators .Ap n (generalizing (2.17)) 



4w = z ± 

and the polynomial 



pz p \ dz 



9 4- 



p— 1 



with[„4± ^] = ±1, (2.49) 



Pkix) :=Y[(x-Wi), P {x):=l. 



(2.50) 



We now have the analogue of Theorems 2^3 and 2J5 , but for the fc- vector p-reduced 
KP-hierarchy: 



Theorem 2.12 The functions ^^ n (x,0; z) in (2. 41) are, for x = 0, the unique 
solutions of the following differential equations 



4=1 j^i 



*± n (0,0;z) = 
(2.51) 

supplemented with the asymptotic conditions (2.48). They are the KP wave and 
adjoint wave functions at t = which correspond, in the Grassmannian settings, 
to the subspaces W^ n defined by 



W± n := span{(^J^±„(0,0;z)}, 



>o- 



(2.52) 



The Wp n satisfies the k-vector p-reduced condition (2.42) in definition 2.10, for 



W p % := P k {Ain)WL C W± 



We also have 



A ± W ± C W ± 

The Lax operators C pn (x, 0), at t = 0, are given by the formula: 



(2.53) 
(2.54) 



4=1 



y± n (x J 0;z) = z p V± n (x,0;z), (2.55) 



and so they satisfy (2.43) and (2.44) w ith 

qf(x,0) = T ^ ±XWl 



(2.56) 
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Proof : At first, check that 



d_ 

2vr J r + dy 




VP+ 1 

e p +1 



-uy 



i=i 




P 

2vr Jv 



e p 




P_ 

2tt /i 



^ y-Wi 
Dul + M + E 



W; 



8=1 







n,; I W; 

cm 



i=i 



8=1 



(£)' +B+ |>( 



9m 



and similarly 



9_V 
dv 



— v 



E 

i=l 



9 

rii | - w, 

ov 



%,ni V )- 



(2.57) 



This yields 



P k ± 



9_ 
"9m 



4Y-. 

aw 



1=1 j^j 



A _ 

"9m 3 



= °> (2-58) 



upon acting on the above two identities with P&(±^-) or alternately one may 



(2.57). 



directly check (2.58) in the course of the above argument and reinterpret it as 



Equations (we recall that D = J|) 



^{±D p ) - x - Y^n i e ±XWi {±D)- 1 e^ xw ^j $± n (x + z p ) = z p $± n {x + z p ) 



follow immediately, upon noting that 

(±D - w)- 1 = e ±xw (±D- 1 )e TXW , 



(2.59) 
(2.60) 



yielding (2.55) and (2.56). 



Also observe that (2.59) is reminiscent of (2.28) and is an example of (2.43), with 



(L +p ) + = D p — x. In the same way as in Theorem 2.6, deduce 
d 



u=zP 



z± n-(^) e±i $^ ((^J^± n (0, 0; z)) , (2.61) 
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and rewrite (2.58) as 



± 



du 



U \Pu ± 



~<9-u 



K ± 



"9m 



E«.n(± 



A _ 

du 3 



(2.62) 

and thus (2.61) and ( 2.62[ ) immediately yield ( 2.51 ). Note that (2.51) uniquely 



define; 



Lemma 



*± n (0,0;z) = l + J^afz-', since (2.51 



2.4 for notation) 



is of the form (see proof of 







z"— + S k -! + ... tf* (0,0;*) 



dz 



and so the corresponding points in the Grassmannian W5 n are uniquely defined 



by (2.51) and (2.52), by almost the same argument as in Lemma 2.4 



We shall now check (2.53); equation (2.51) implies 

^(4)t(0,0;z)eW± 



(2.63) 



while for t > 0, 

z p P k (A± n )(A£J*± n (0, 0; z) = z p (J^JP k (^ n )^ n (0, 0; z) 

= V + (^n)']) ^(^n)^n(O, 0; *) 

MJJ^ P + E (^n)i^^^J(^n) J V(^n)^n(0,0;z) 
i+j=£-l / 

= ((^n) V - ^J- 1 ) ft(^J*J„(0, 0; *) G W p ± n , 

from (2.49) and (2.52) upon using ( |2.63 ), yielding (2.53); also the codimension 
statement is obvious as well as (2.54). To see that ty^ n (x, 0; z), as given by (2.47), 



are indeed the wave functions at t = 0, by precisely the argument of Lemma 2.5 
we need to show (see (2.22)) 

^Wtf 1 - ^) = ±D*± n (x,0;z), 

which follows from (2.47), using (2.46) and ( 2.61[ ), concluding the proof of the 
theorem. ■ 



6 Using 5 k -i{l/z e ) ^ z k £(l/z l ),£> 1 and S k -i{l) = 0. Note in Lemma 2.4 k = 0,5-! = 0. 
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3 p— Airy kernels, vertex operators and Virasoro 



In this section we write the basic kernels of our theory in terms of KP wave opera- 



tors, leading to Virasoro identities for the associated Fredholm determinants (3.9) 



and various useful identities (3.7), (3.16) and (3.48) for the kernels, expressing 



them in terms of wave functions or double contours integrals. 
Vertex operators typically generate Darboux transformations in integrable sys- 
tems at the level of the tau functions. The KP vertex operator X(t, y, z) is defined 



as: 



X( t ,y,z) := -J—e^-y^e^"-^ 
z-y 



1 JL 

at,. 



(3.1) 



Given two integers p, n > and w G K let us consider the KP wave functions at 
t = 



*p, n (x,0;z) 



±p it 



2tt 



-Z 2 



,P+1 



(y-w) ±n dy, (3.2) 



where the paths are chosen appropriately as described in the previous section. 
According to sub-sections 2.2 (for the case n = 0) and 2.3 (for the case n > 0) 



the objects in (3.2) are wave functions for the p-reduced KP hierarchy and for 



the 1-vector p-reduced KP hierarchy respectively. The corresponding (pseudo)- 
differential Lax operators, at t = 0, read 



C*(x, 0) := (±D) P -x- n(±D - w 



,-i 



[±D) P -x- ne ±xw (±D)- 1 e Txw . 



(3.3) 



In the sequel, when the subscript n is omitted, we mean n = 0, coherently with 
the notation in the subsections 2.1 and 2.2. 

We can now define the basic objects of our theory, remembering D^ 1 means J dx. 



Definition 3.1 Given the wave functions (3.2) we define the integral kernels 



KK(A,A') 



e p +1 



p_ z p+i k x ^(Z)Z >s ) v z ip+i 



2TTpZ P -^ +n z' P -^- n 



QP+l 



(3.4) 
(3.5) 



*=aVj> 



Also let E := Ui=i[ a 2i-i) Q2i] C M + be a disjoint union of intervals and let 

Ep := {x E R + such that x p G E}. 
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We define, for a tau function r(t) ; another function T^(t) 

-\x \ dz X(t; ojz, J z) 



r E (t):=e Je1/p r(t). (3.6) 



The function Te(t) will be shown to be a tau function itself in Proposition |3.3| 
The following propositions, due to Adler, Shiota and van Moerbeke in jl], will be 
crucial. In the following we denote (: : means normal ordering) 

WP = A + (-i)t-i, } = £ : wf >wf > : -(t + > 

1 i+j=e 

o p 2 — 1 

and cf. ] = 5\ » r-. This first proposition is a consequence of Theorem 6.1 of 

I2p j 

HI- 

Proposition 3.2 T/ie kernel k~\(z,z') satisfies the equation (see (2.4) for the 
notation i ) 

r(t) 

where we have an identity of formal expansions in Sato's theory. 

This second proposition explain how to express the Fredholm determinant of 
K^l in terms of the KP tau functions we described here and in the subsequent 
sections. 

Proposition 3.3 Consider a disjoint union of intervals E := IJi=i[ a 2i-i> a 2i] C 
R + and the Fredholm determinant det(l—27r/iK^'lx E ). Given oj ^ u' two p-roots 
of the unity, the following equality is satisfied 

det(M-2n^X E ) = ^. (3.8) 

Moreover te, as t, is a k-vector p-reduced KP tau function (k = 0, 1 when 
n = 0,n > respectively) and the following Virasoro constraints are satisfied: 



r(t) 

T E (t) 



2r 

a 



(3.9) 
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where c^j are constants, given above for the case n = 0. These constraints are 
valid, for n = 0, Vj > while, for n > 0, they hold just for j = 0, 1. 

Sketch of the Proof: 



Formula (3.8) is just Corollary 7.2.2 of [4J, modulo the conjugation term 



e p+i z " z ' n in K^fl , which has no effect on the Fredholm determinant. 



L X,t 1 



z~ e p +1 ' 

Now we have to prove the Virasoro constraints and the fact that te is a tau 
function. 

First let us consider the n = case. In order to deduce (3.9) from the Corollary 
7.2.2 of [I], the first fact is that by the invariance conditions (2.18), we find that 
the point in the Grassmannian Wp satisfies z^ p ApWp C W+, with 



z 3P At 



y (j-i) p +i_9_ 



P 



+ z 



jp+1 



p~ 1 Jj-l)p 



dz ' 2p 
and so by (2.15), one has (V + f .+)- = 0, and in particular 

1 



j>0, 







P 



(l+) 0-i)p 



p • 



which by the Adler-Shiota-van Moerbeke correspondence [3] leads to 



1 



2p 



W$L l)p + ($j T (t) = 0, forj>0. (3. 



10) 



(Here and in the sequel we denote with the subscripts "p" or "p, n" the Lax 
operators related to the particular solutions of the relevant KP hierarchy we are 
dealing with.) The first term in the operator contains ^{p + l)t p+ id/dtj p+ i and 

thus the shift t p+ i i-> t p+1 — has the virtue to eliminate the W^ +1 term in 
(3.10), thus yielding 



( L w w _ p —w [l) + c (0) i t o 



Both, Corollary 3.2.1 and Theorem 4.1 in [4], and u p = u' p = 1, yields 



d_ 

dz 



[z^-^ p+1 X{t,uz,u'z)) 



1. 



(2) 



P 



-W, 



(i) 

0'-i)p 



i (0) 



X(t, coz, u'z) 



without the W^-^ -term.. It is legitimate to add this term, because 



(3- 



(j-i)p 



(v* - u e )X(t, u, v) = [wf \ X(t, u, u)] 
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and thus 



[wj 1 ^, X(t, u, v)] = for u = uz, v = u'z andp|£. 



Then doing the shift again t p+1 i— > t p+1 + reintroduces the Wj^-term again 
in (3.11). From the identity (3.11) and (3.10), it then follows from the arguments 



of Theorem 4.1 of [3] that te defined by (3.6) satisfies 

2r 







u- p - 1 w m 



c^)\r E = 0. (3.12) 



2p ' ' U-i)p P>i 

That te is actually a r-function follows from several important facts 

(i) if r is a r-function, then e aX r = (1 + aX)r is as well. 

(ii) the vertex operators for different indices commute: 

[X(t, A, /x), X(t, it, t>)] = 0, for u 7^ /i, A ^ v. 



(Hi) The integral in the exponential (3.6) is the limit of a Riemann sum, which 
using the higher Fay identities guarantees that t e is a r-function and it is ex- 
pressible as a Fredholm determinant, as explained in [3]. Moreover, for n = 0, r 
and te are p-reduced KP r-functions; indeed r is, because the above satisfies 
z p Wp C W+. Also r E (t) is, because X(t; uz, w'z) is missing t ip and J^-, since 



1 for z > 1, concluding the proof for the n = case. 



Now for k = 1 (i.e. n > 0), we have that the 1-vector p-constrained tau func- 
tion is also characterized by the fact that r and (d/dt p )r are both tau functions 
j. Since the vertex operator in (3.8) is free of t ip , (d/dt ip ), i G Z + , (3.8) yields 



E 1 /p 



dz X(t, wz, w'z) q t 



dt p dtp 

and this is a r-function since (d/dt p )r is a r-function and e" M ^ 1/pdz x(t,w«X«) 
takes r-functions to r-functions, as mentioned before. In particular, therefore te 
is a r-function as well as (d/dt p )rE, and hence te is a 1-vector p-constrained 
r-function. 



To see (3.9) for the n > case we follow roughly the n = argument. By (2.42) 



(2.53) and (2.54) (we recall that here k = 1), one has the following inclusions: 

1 



A+ n W+ n c W+ n , z*(A+ n -w)W+ n c W+ n , for^+ n (^) 



d p — 1 
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By (|2.15b, one has (V 



0, and thus one has (3.10) for j — with 
p — 1 



i — y n — 



r p — 1 



2 v 2 7 

This proves, by the same argument as before, the identity (3.9) for j = 0. 
For the case 7 = 1, the relation follows from (P, + „. .)_ = 0, and so we have 
one extra term, not present in the k = case, due to the presence of wz p , yielding 
(3.9) with k = n = for j = 1, with one additional term, 

2r 



£ a -k + k K ' + ~ wW ' 1} + ) T=0 ' 



(3.13) 



It suffices to show (d/dt p )r = — (d/dw)r to conclude the proof of (3.9). To see 
this, notice that the integral ^p n (x, 0, z), as in (3.2), readily satisfies 



dw 



(x, 0; z) = -n(D - w^V+Jx, 0; z) = (L+) p _*+(x, 0; z), 



and upon differentiation by x, 



d fd 



y*+ n {x,0;z) = (L^)_(^)^+ n (x,0;^), 



dw \dx 

and so, acting on W+ n , as defined in (2.52) 

dw p ' n 



0. 



Now let's recall the standard Sato formula (2.9), here rewritten as 



* ± (x,t:z) = e ± ( xz+ ^ t * zi y- 



r(t) 



with r] := 



z 1 d 
i dt 



(3.14) 



(3.15) 



The formula (3.14) leads to 





(e-"-r 



dw 



T 

V / 



dw 



^p,ni x ^ z ) (L+p) fact: 

y ' V p,nJ— p,n\ > ) 



(e^-1) 



at. 



r 

V / 

using in the first equality straight differentiation and in the last equality again 
the Adler-Shiota-van Moerbeke correspondence j3], (remembering rj from (3.15)). 
This shows 

dr dr 
dt n dw ' 



as claimed, concluding the proof of the Proposition 3.3 
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Proposition 3.3 then leads to the following theorem involving bona fide ana- 
lytic identities, rather than identities involving just formal series. In the following 
we denote with the subscript t> p= o the locus given by {U = Vz > p}. 



Theorem 3.4 Setting all t{ = 0, Vz > p, one has the following kernel identity 

e -V p (u)+(z'P+x)u / u _ w - 



for the kernel defined in (3.5), 



i 



du I dv- 



-V p (v) + (zP+x)v \ v 



where V p (u) is a polynomial 



w I u — v 

(3.16) 



V p (u) := — — + V 6— — , ^ := 0^, . . . , 
p + 1 r-rf « + 1 



(3.17) 



i=0 



wzi/i 6*j polynomials implicitly given in terms ofti, ...,t p _i, by solving the equation 
w = V'(u) for u in terms of a series in large w, as in (i) and identifying it with 
another series, as in (ii); thus 



u 



(0 



(n) 



ujp l_ 2 w p dp-Z w p 

p p p 



1 / p-3 



0l_ 2 + 6> 4 w p 



2p ~ p - 2 " 
H h 0{w^p) 



(3.18) 



w 



1 

» + - ^(p - 3%-jW-i + 0(w- 1 - 1 / p ). 
P i 



Example: The first 0j satisfy the following equations: 

1 



p — 1 



-ft 



P-2 



p-2 



-0 



1 



p-3 



p — 3 



-0 



p— 4 



The proof of Theorem 3.4 requires Proposition 3.6 which itself is based on 
the following Lemma: 

Lemma 3.5 Consider the solution to the p-reduced KP-hierarchy coming from 
Wp C Gr as in Theorem 2.3. The associated differential Lax operator Cp~(x,t) = 
Lp P (x,t) satisfies 



c;,(i + p -±^t p+1 ^D+ L - J2 kt k (c;' p 



k—p 



fc=2+p 
k^ip 



and thus for t> p = 0, one has 

[c;{x,h,...,t p ^),D] = i. 



(3.19) 



(3.20) 
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Proof of Lemma 3.5: For notational convenience let's denote 



(3.21) 



with W p the dressing operator as in (2.6) associated to this particular solution. 
Our first claim is that the corresponding operator M p of (2.7) satisfies the equa- 
tion 



k=i 



Indeed, from (2.7), one has 



M; = W p xW p ~ l = S p e^^ Dk xe-^^ Dk S;\ 
Since [D,x] = 1, one has [f(D),x] = f'(D), and thus 




(3.22) 



(3.23) 



dU 



and so 



+ Y,kt k D k -\ 



(3.24) 



k=l 



Setting this formula into (13.23) and using S P D 3 S p 1 = (L+ ) J yields (3.22 



The conditions z p W p C VV p and A p W p C implies, upon using (2. 



5), that 



both Cp = (L+) p and are differential operators. One then computes, using 



fl2.14| |, ( |2.15D and (J3_22j) in the third line, 
1 = [.4+ = = [(P+), (P+ ± ) + 



2p v p 



£+ (S p DS p 1 ) + + -(S.xD^S; 1 ^ + - J>* {( L D k ~ P ) + - ^(S P D- p S p 
p p k=i ->' 



fc=p+l 



\ " / " fc=p+2 



(3.25) 



yielding (3.19) and hence (3.20), upon setting t>„ = 0. 
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Proposition 3.6 Consider the Lax operator £+ (x,t) flowing off the initial con- 
dition Cp (x, 0) in (3.3) together with its wave functions \l/^ n (x,t; z) flowing off 
the initial conditions W^ n (x, 0; z) given in (3.2). Then, evaluated at t> p = 0, 
these quantities have the following form: 



(3.26) 



where V p (y) is defined in (3.11), and 



2tt 



^v p (y)±(x+znv( y _ w )±"dy, 



where the polynomial P(t) = P(t\, . . . , i p -i) is determined by 

d -P(t) = {(v;(y)) i/p )_ 



(3.27) 



dU 



1 < i < p - 1 and P(0) = 0. 



Proof of Proposition 3.6: 



We divide the proof in three steps; in the first step we prove (3.26), for the case 



n = 0, in the second step we prove (3.27) again for n = while in the third step 



we extend both formulas to the case n > 0. 

Step 1: Let n = 0. Since Cp(x, ti, i p -i) is a differential operator of order p 
satisfying the commutation relation (3.20), namely [££(x,ti, ...,£ p _i),D] = 1, it 
must have the following form: 



p-2 



fc=0 



£ p + (x,ti,..,Vi) = (L+ p ) + = J D p + ^^(ti,-,Vi) J D fc -x=:g( J D)-x, (3.28) 
and thusQ 

((^r /? ) + = ^+E c «(^-'Vi)^=((Q^)) n/p ) , i<n<p-i, 

i=0 + 

(3.29) 

with the c ni being polynomials in the 9j. One concludes that 
dQ(D) dC+ 



n-2 



dt n 



dt r 



m;r /p h,^} = kq(d)*) + ,q(d)-x] 

= [(Q(D)p) + , —x) 

d(Q(D) n '*)+ 
dD 



7 Since (£+) p = L> + J]^ 2 Oi(t)ZH + a p (t, x)D-p +1 +...,(£$ ) + is free of a; for n < p. 
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and thus, replacing D by y, 



dQ{y) _ d 
dt n dy 

Next, for large y, solve the equation 

P-2 



{Q(v)>] 



Q(y) = y p + J2 e ^-> t p-i)y k 



k=0 



implicitly for y = y(Q,t), namely 

i 1 _i 1 _2 1 

y = Q p — Op-zQ p — P -3Q p — 
p p p 



dp- 4 — 



PZ3 
2p 



p-2 



Q~ p 



Then, using (3.30), one computes 

dy dQ dy 







dy 
dt„ 



dQ dt n 9t n 

n 

dy d(Q(y) p ) + 
dQ dy 



+ 



dy 
dt n 



dy dQ(y)p dy_ Q (—(- 
dQ dy dQ \dy \y 

— / 



+ 



dy 

dtn 



dQ dt r 



dy 



-A-l 



and thus 



^ = !^g ! T £ + 0(g- 1 -?), l<n<p-l; 



<9f 



P 



hence, y as a function of Q has the following form: 

p-i 



(3.30) 



(3.31) 



+ 0(Q- 1 -?). 
(3.32) 



(3.33) 



(3.34) 



Equating (3.32) and (3.34), we solve for 6^ = 0j(ti, 1 < z < p — 1 
inductively as polynomials in the t, yielding = Vp(y), where V p (y) is defined 
in (3.17) and (3.18) of Theorem 3.4, concluding the proof of expression (3.26) in 
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Proposition |3. 6 



Step 2: To show that ^ = \l/ p0 defined in (3.27) are the KP wave functions, it 
suffices to show that they satisfy the first equations in (2.8) and in (2.10). Let us 
first consider the former for ^ p . Indeed, the integral 



F P (X) :- 



(3.35) 



satisfies 



d_ 
r+ dy 



V p (y)+Xy dy= / (_y'( y ) +X ) e -V p (y)+>v dy 



-V' 



d_ 
OX 



+ X)FJX). 



Using the form (3.26) of C p (x, tx, i p -i) with n = 0, it shows that, setting 
A = x + z p , 

C;(x,t u t p ^)F p (x + z>) = (V;(D) - x)F p (x + z*) = z p F p (x + z p ), (3.36) 

establishing C p (x, ti, t p -i) ^ p = z p ^ p for ty p as in (3.27) with n = 0, and 
similarly for establishing (2.8). 

We now establish (2.10); from (3.35) and using (3.30), namely J^(gf~^p(y) + 
{Vp(y) n ^ p )+) = 0' integrated from to y and remembering V p (0) = 0, compute 



d_ 

dt r , 



F p (x + z* 



V p (y) ) e - v ^ y)+(x+zP)y dy 



d_ 

dt n 



y=0 



e -v P (y)+(x+z p )yri y (3.37) 



(v;(D) a P ) + F p (x + z*) - (ij(y)5; 



y=o 



F p {x + z p ). 



The expressions 



h{t) := {V'{y)i) + 



y=0 



1 < % < p 



(3.38) 



j tp— 1 ) 



are polynomials in ti 

by the KP flow compatibility conditions and by (3.29) that for 1 < i, j ] <p — 1 

d 



dt 3 



as follows from the proof of Lemma |3.5| Observe 
lit ions ai 







dt 



{£>] 



d_ 
dt. 



(4: 



(K(D)r 



\D=0 



l_D=0' 



(££)+,(£ 



lD=0 



(K(D)i) + ,(v;(D)i] 



\D=0 



0. 
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Thus, the polynomials bi(t) can then be represented as 6j(t) = J^P(t) for some 
uniquely defined polynomial P(t), modulo an additive constant, which we set to 
= 0; so, we have 



be 0, i.e., P(0) = 0; so, we have that expression (3.38) equals 





=o dt 



P(t), l<i<p-l. 



(3.39) 



Setting this expression into (3.37), and defining ^p(x, t±, £ p _i; z) as the right 
hand side of the first expression (3.27), leads to the following differential equations 
satisfied by ^j" for 1 < n < p — 1: 



dt n 



^-z 2 e p+i z —(e p ( t >F p (x + z p )) 

Z7T Ot n 

P 2=1 £—zP +1 

-z ■ e P+1 




^ e - 5 fr^ +1 (£+p) + ( e ^)F p (x + ^)), using gg} 

(3.40) 




P sz 

2^ 



with the initial condition, by (3.2) with n — 0, 



*J(s,0;z) = A /f ^e-^ 2P+1 e^F p (x + 

V 27T 



*")|t=o 



thus establishing (2.10) for and similarly for , and thus we find the identity 
(3.27) for n = and conclude the proof of Proposition 3.6 for n — 0. 
Sfep 5: The proof for the case n > follows from that of n = with some 
important modifications. Indeed, since V% = IZ n ^ (L£ n ) + , we find as in 
Lemma 13.51 that 



1; •••) vp—l)-i 



D] 



but with 



L£ n (ti, V-i) = Vp{D) -x-n{D- 



w 



(3.41) 



(3.42) 



with V p (y) specified in Theorem 3.4 Equations (3.41) and (3.42) are proven 
exactly as before. Indeed, one sets, analogous to the n = case, with Q as in 
<p8b, 



2£ B = Q(D) -x + (2£J_, and so, (Z> n ) + = (Q{D)\ 



l<j<p-l, 



and the same proof as before yields Q(D) = V'(D). Indeed, one observes that 
(<2(£>) i/p )+] = 0, 1 < j < p- 1, since (Q(L>) J >) + is x-independent; thus 
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we have: 



dU 



L p p , n ,(Q(D)- P ) + 
(L p P , n ) + ,(Q(D)t) 



+ 



showing that d(L^ n )^/dtj = 0; so ( L^ n )_ is unmoved by the first p — 1 flows 
(but moved by the t p flow), yielding (3.30) as before, etc. and so formula (3.42) 



is indeed proven exactly as in the case n = 0, while (3.27) is proven exactly as 



for the case n — 0, completing the proof of Proposition | 3.6 
Now we can prove Theorem |3.4| 



Proof of Theorem 3.4' First consider th e ca se n = 0. Setting A(z,z') : 



2irp z'^ep+ lZ e~p+ lZ '" and using (|3.7|) together with the definitions in 
D (a(z,z')K^(z p ,z'^ 



(PI 



(3.43) 



M — V 



Then using J^e 11 ^ ^ 
and (3.43) above, one checks 

A(z,z') 



]e x(u v) ^ representation (3.27) for ^ = \1> 



± 

p,0 



(2vri) 



du / (if 



; -V p (u)+(z' p +:c)« ^ 



1 P 



A(*,*0 
(2vri) 2 



-Vp(u)+(zP+3:)u 
e V 9 (v)-(x+*>)v dv 



-V p (u) + (x+z'P)u du 



r+ 
1 p 



*> P =o 



*> P =o 



using (3.27), 



using (3.43), 



thus yielding the desired identity (3.16) of Theorem 3.4, except for the differen- 
tiation D. To do the identification, without the D, one hits these identities with 
D~ l , sets D~ l e x<yU ~ v ^ = — and then identifies the formal expansions (of the 
Sato theory) of the wave functions ^ of Proposition 3.3 (which can ultimately 
be traced to Lemma 2.5) with the asymptotic expansions of the double integral 
(which can be ultimately be traced to Theorem 2.8), which is an analytic ob- 
ject, concluding the proof of Theorem 3.4 for n = 0. This proof can be adapted 
without trouble to n > 0. ■ 



Example: For p = 2, 3 we find from Theorem 3.4 that 
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u 3 

V 2 {u) = — -t\u 

Vz(u) — — — t 2 u 2 - t x u. 
Hence, upon choosing the appropriate contours as in Corollary (2.9), we rec- 



ognize, for n = 0, the well known expressions of of the Airy and the Pearcey 
kernel as double contour integrals. The case n > corresponds, for p = 2, to the 
case of the kernel for the Airy process with outliers, studied in [5] and [6j. For 
p = 3 we recover the kernel for the Pearcey process with inliers, again studied in 

m- 

Lemma 3.7 For n = the r -functions corresponding to the points in the Grass- 
mannian W + such that: 

evaluated on the locus ti — 0, i > p, i.e. the so-called topological tau function 
Ty P \ti, t 2 , tp-i), are completely determined alternatively 

— — logr (p) = res u (V;(u))*, 1 < % < p - 1, (3.44) 

and the KP hierarchy equations containing just di, d 2 , ■ ■ ■ , d p , or equivalently 

— logr (p) = P^ Te s u (V;(u)) P -^ , 1 < % < p - 1, (3.45) 



with V p (u) given in Theorem 3.4 So in particular we find: 



logr (2) (t 1 ) = -^?, 
logr (3) (t 1 ,t 2 ) = -^t 2 -^4 (3.46) 



We now give another representation of the kernel in terms of the wave and dual 
wave functions: 



3 The residue is evaluated at u = oo, i.e. it is the coefficient of u 1 , for u 1 small. 
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Theorem 3.8 Setting all £j — OA > p, one has the following kernel identities 
for the kernel defined in (3.5): 



t> p =o ~~ z'p - zp 



\k+e=p-i 



i=l fc + £=j_i / 



i> P =0 

(3.47) 



T/ie pref actor in (3.4.7) reads 

<p(z,z?) 



■ J T (2'P+ 1 -zP+ 1 ) / ,\n 

ie^ +1 i z ^ 



2irp(zz') P, 2 V : 



and ^p n (z) = typ n (x, t, z)\ t>p =o are given explicitly in Proposition 3.6 (the de 



pendence of^^ n on the variables (x,t) has been suppressed for notational's 
venience). 



con- 



Theorem 3.8 is an immediate consequence of (3.27) and the following useful 
lemma. 



Lemma 3.9 Let us denote 



2ni / r ± 



Then we have 



v ~'p\ 



p-2 



1 



Z'P - ZP 



p^+mx-^v, 



(3.48) 



yk+t=p-l 



+ E {D k ^{z^-DYr p {^+n{D-w)-^{z'){-D-w)-'r p {z) 



i=l k+t=i-\ 
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Proof of Lemma 
1 





3.9 



Observe that, because of our choice of F^, 



(2m) 
1 



, / ■ > 9 9 

du I dv I - — h 



du dv J V v — w 



u — w 



-V p (v)+{x+z'*)u i 
-V p (v)+(x+zP)v u _ v 



, , du I dv 

(2m) 2 U J r - 



-V p (u)+(x+z'P)u 



-V p (v)+(x+zP)v \ v — 



U — W 



W 



■Vp(u) + v;(v) + z'p - z p 



+ 



(2m) 2 
z f P _ zP 
(2m) 2 



du I dv- 



-V p (u) + (x+z'P)u 



r- c 
1 p 



-V p (v)+(x+zP)v \ v — 



U — V 



u — w 



n 



(u — w)(v — w) 



du dv- 



-V p (u)+(x+z'P)u 



1 p 



w 



u — w 



v;(u) - v;(v) 



u — V 



+ 



n 



(u — w)(v — w) 



Vp(v)+(x+zP)v \ v _ w l u — v 



1 



The proof is finished upon using (3.16) in the last expression and noticing that 

e -V p (u) + {x+z'P)u / u _ w 



by 03.17) 
1 



(2m) 



du / dv 



p " *- p 



(2m) 2 
1 

2tH 



du / dv 



Vp(v)+(x+zP)v \ v _ w 

e -Vp(u)+(x+z'P)u / u _ w 



V'(u) - V'(v) 



V p (v)+{x+zP)v \ v — 
1 



u — V 



+ 



p-2 



n 



(u — w)(v — w) 



u — v 



+ 



n 



du e-^(")+(^ z ' p )> - w ) n 



ir I \ a — r * — ' u — v (u — w)(v — w) 

i=l v ' v ' , 

e Vp(v)-(x+zP)v 



2m 



dv- 



[v — w 

p-2 



i=i fc_)_£ = j_i 



(u — — u>) 



p-2 



i=l k+£=i-l 



Example: Let us set n = 0; for p = 2, 3, one uses again the contours chosen 



in Corollary \2.9\ we recover using Lemma |3.9| the expressions of the Airy and 
the Pearcey kernel as integrable kernels a la Its-Izergin-Korepin-Slavnov [24J. 
Indeed, for p — 2, denoting with Ai(x) the usual Airy function, we find, setting 
£>/ = /', that 

.( 2 ) M xrv _ Ai'(x + ti + X')Ai(x + t 1 + A) - Ai(x + *i + \')M'(x + ti + A) 



A' -A 
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which is precisely the Airy kernel for x + 1\ — 0. For p = 3 let us define 

which is nothing but the so-called Pearcey function together with its adjoint. We 
obtain 

41* 2 (A,A') = 

ij)$"(x + h + A', t 2 )^3 (x + ti + A, t 2 ) - ipt'( x + ti + A', t*)ij)z{x + h + A, t 2 ) 

+ + ti + A', ^^"O + *i + A, * 2 ) - 2t 2 ^(x + ti + A', t 2 )^ 3 ~(x + ti + A, t 2 ) ) (A' 

the usual Pearcey kernel with parameter t — t%, for x + 1\ — 0. 

4 PDEs for random matrix kernels 

The machinery of the previous sections will be used to compute PDEs for Fred- 
holm determinants of kernels occurring in random matrix theory and Dyson Brow- 



nian motion, by specializing the kernels in (3.5), thus yielding a proof of Theorem 



The PDEs will come from PDEs of the KP hierarchy exploiting the formula 



( ]3.8[ ). Then, using the Virasoro relations (3.9), we will replace some of the partial 
derivatives in ti with partial derivatives in a%, i.e. the end points of the disjoint 
union of intervals E := UI=i[ a 2i-i! a 2i] C 1R. The next lemma provides the PDEs 
of the KP hierarchy given in the Hirota form which we shall use to generate the 
promised PDEs for the Fredholm determinants; we will denote with pi the Schur 



polynomials already defined in (2.5). The following lemma can be deduced from 
arguments in [20] . 



Lemma 4.1 The bilinear identity for KP (2.1) generates two strings of Hirota 
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relations^ 

dz 



} T{t-[z-^)r{t' + [z^])e^^^\^ 
2m 1 ,, 



= Er=oPi(^+i(^)e^^^Tor 
= IX i Jfc {Pe+i(dt) - l<hdi) tot 

+ E^ 2 yivt-i (pi+i{dt) - \d2dt-1 - \dm{d t )) ro T + o{yf), 

which are independent, for I > 5. The first string, denoted symbolically by Y^ ; is 
the standard KP hierarchy and we will denote twice the second one minus twice 
the first one by Y lj ^_ 1 

Y e : (p i+ i(d t ) - -dxd^jro r = 0, Y^-i : (didi> - -d%di-x - 8^(8^ r o r = 0. 

(4.2) 

We are ready to state the main theorem of this section; we remember that, given a 
disjoint union of interval E := Ui=i[ a 2i-i> a 2i] C K, we defined in the introduction 
the two differential operators 

i i 

Also we will denote 

Q = Q P (h, ...,t p ^;E) : = logdet(/ - K^X E )\ x=u=t>p=Q . (4.3) 



Theorem 4.2 Each of the Hirota equations 

Y 3 , . . . , Y p+1 , Y M , . . . , Y 1|P and 2(p + 2)Y p+2 + (p + 1)Y 1)P+1 , (4.4) 

gives rise to a non-linear PDE for the Fredholm determinant Q. These PDE's 
only involve the differentials 8 and e with regard to the boundary points of E and 
the t-partials 82,83, <9 p _i. 

9 We recall the standard notation for the Hirota symbol of two functions / and g, associated 
with any polynomial of many variables 

( d d \ 

p{d u d 2 , ■ ■ .)fog :=p I — , — , . . . I f(t 1 +y 1 ,t 2 +y2, ■ ■ -)g{h - Vx, h - 2/2, ■ ■ -)| {l/i}= o 
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Examples of PDEs are given below; also note that linear combinations of these 
equations can lead to new interesting equations. 

Corollary 4.3 Q satisfies the following PDEs: 



Case 1 (n = 0) : 



• For all p > 2, Q satisfies 
the Y 3 -equation: 

d 4 Q+6(d 2 Q) 2 + 5 2 , p (2-4e)dQ 



+ (1 - S 2>p ) 3d 2 2 Q - 4 



£—^tr-id + (l- 5 3jP )d 3 
p 



• For all p > 3, Q satisfies the equations (4-5) and 
the Y 4 - equation: 



(4.5) 



d 2 d 3 Q + 6(<9 2 <9Q) (d 2 Q - -(p - ) + - Se)dQ + 2t 2 dd 2 



p 

+(1 - 5 3}P ) \2d 2 d 3 Q - 3 
the d 2 Y 3 — dY \- equation 



V 



(p-2)t^ 2 d + (l-5 4;P )d 4 



did 



o- 



8 - 2(1 - 6 3 , p )d 2 d 3 dQ + 5 3tP (e - 2t 2 d 2 - 
- -(1 - 8 3tP ) ((p - l)t p ^d 2 - 2(p - 2)t p _ 2 d - |(1 - 5 AiP )d} d 2 Q = 0, 

the d 2 Y 3 -equation (Boussinesq form ofY 3 equation) 



d A U + 3d z 2 U + 6d 2 U 2 - 4(1 - 5 3iP )d 3 dU = 0, U := <9 2 Q - 

• For a// p > 4 ; Q = Q p satisfies all previous equations and 
the Y 5 — Y lj4 - equation: 

d 2 d 2 Q + ld 3 d 3 Q + ^ 3 2 Q + A(d 2 Q)(d 3 dQ) + A(d 2 dQ) 2 



p — 1 
p 



i p -\. 



+ 2 is* 

16 

p 



P 



p+2 



4 



p — 3 



(4.6) 



(4.7) 



(4- 



d 2 res u (V'(u))— (d 2 Q) - -(p - 1)<9 3 <9Q - 12^ *p_ 3 (l - K 



p + 2 ' r ' J p 

(p - 2)t p _ 2 d 2 dQ + <5 4 , p [1 - 4e + 2t 2 «9 2 + 3t 3 ^.! 



P 

+ 4(1 - £ 4 ,p)(l - S 5jP ju 5 
(4.9) 
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Case 2 (n > 0) 



p = 2 : <9 4 Q + 6(<9 2 



(2 - 4(e - wd v 



+ Sd 2 



(4.10) 



p = 3 : d A Q + Q(d 2 Q) 2 -8t 2 d 2 Q + 3dlQ + 4d v 



(4.11) 



(<9 2 <9 2 - 2t 2 <9 2 - 3(e - wd w ) + l)dQ + 6{d' 



2d,d, 



(4.12) 



(e - wd w - 2t 2 d 2 - 2)d 2 Q + dlQ + 2{«9 2 <9Q, d 2 Q} 9 + 2d 2 d v 



(4.13) 



<9 4 £7 + 3d 2 2 U + 6d 2 U 2 + 4d w dU = 0, 



U = d 2 



(4.14) 



Notice that for n = 0, the function Q satisfies these equations, but without the 
terms containing d w . 



Note that Theorem |4.2| and Corollary 4.3 immediately imply Theorem 1.1 
Proof of Theorem \4 . E " We first give the proof for n = 0. Notice that 



TOT 



0, for J^ii odd. 

i 

j J log r + other partials of log r, for = even. 

i i 

(4.15) 



Taking into account this remark, one has for the Hirota equations (4.4): 
(i) The partials d p can be ignored, since the r-functions appearing in 



det(I-K<*lx E 



r E (t) 

T(t) 



are r-functions for the p-reduced KP-hierarchy; that is the r-functions do not 
contain t p , t 2p , . . . and thus d p = 0, when acting on those functions, 
(ii) The odd-degree terms in the Hirota operator (4.4) do not matter, because 
they vanish as an Hirota symbol acting on r o r, by (4.15). 



(iii) The Hirota equations Y p+1 and Y l p contain d p+2 and did p+ % and many other 
terms involving di for 1 < i < p — 1. The partial d p+2 can be omitted, since it is 
of odd degree, and d\d p+ \ will be taken care of. 

(iv) The Hirota equations Y p+2 , Yi )P+ i contain d p+ 3,d 2 d p+ i, did p+2 and d 2 d p+ \. 
The first two terms do not matter, since they are of odd degree, and the exact 
linear combination 2{p + l)Y p+2 +pY 1)P+ i removes the term did p+2 , but d 2 d p+ i 
will need to be taken care of. 
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To conclude, the list of Hirota equations (4.4), acting on p-reduced r-functions 
only involves di, . . . , <9 p _i and derivatives with regard to higher ij's only through 
did p+ i and d 2 d p+1 . 

From the Proposition 3.3, in particular (3.9 ) for n = 0, both g : = log T^(t), logr(t) 
satisfy 



dg 



£9 



f - ^ itj + <9i j # + T p =: dg + T p 

V i>p+l / 

ili di + d P+^j 9 + c p =: Sg + c p , 



where (refering to ( |3.9[ )) 



(0) _ P 2 -l 
12p 2 



and of course <91nr = 0, elnr = 0. Using [di,5] — for 2 < i < p — 1, compute 
inductively 

<9 m s = c>% = d l (5g + F p ) = 5d*g = 5(5 l g + S^%) = 5* +l g + 8%, 

n ^ = n + ^-%) = & n + ^ n 

i=2 i=2 i=2 i=l 

£^ = e(8g + T p ) = cfeg = S(5g + c p ) = 55g, 
d 2 eg = d 2 (5g + c p ) = d 2 Sg, 
and thus on the locus £ := {t\ = 0, t p = t p+ i = ... = 0}, 

= <9i<9 p+ i + ~ I 9 i + /jtid i d 1 J ,dxg = dg - T p \ c . 
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So we conclude from the formulas above that, on C: 



p-i 



p-i 



Ufa = nafaS-^^PE^^) 

P V 2 i+j=p I 



d 2 d p+ ig 
dig 

d\g 
d\g 



1=2 



*i=0 



(4.16) 



p-i 1 
pdid p+ ig = (pe - l)dg - ^ iUd^g + — {iU){jtj) 



2 

P-I 



i 1 
(pe - -^itiididg - -(p - z%_j) + ^ ( l '*<)0'*i)> 



P 



1 / \ 

<9 2 £# 2<3 2 + y^iUdA g. 

P\ 2 J 



9 2 5- - (1 - <S 2 . P )- — -tp-i, 
V 

d 3 g- 6 -f, d\g = d l g, i>A. 



(4.17) 

(4.18) 
(4.19) 

(4.20) 
(4.21) 



Substituting (4.16 )-(4.21 ) in the explicit PDEs Y e or Y M _i in g = \nr E (t) yields 



an explicit PDE in 5 2 , ds, . . . , dp-i, d, e, for both logr^t) and r(t); since 
logr^ — logr, one finds two PDEs: 

r (log r(t)) = 0, r s (logr B (t))=r E (Q + logr(t)) = 0. 

Then form the PDE 

T E (\ogr E (t)) - r (logr(t)) = r £ (Q + logr(t)) - r (logr(t)) = 



p-i 



in which JJc^* logr(t) remains; then use (3.45) to explicitly substitute its value 



as a polynomial in £ 2 , tp-x- Carrying out the program in a few cases yields the 



equations in the corollary 4.3, case n = 



For n > 0, the only change in the Virasoro, in comparing (3.9) for n = and 
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n > 0, is that r p i— > T p + nt p + c^g, and e — >£ — w-J^ = e' and c^j will have 
a different value, and these changes have no actual effect beyond changing e to 
e'; then carrying out the program will add some extra terms to the equations. 
This proves the PDEs for E C M + , and then we extend the PDEs to E C R by 
analytic continuation. ■ 

In the following appendix we go into more explicit detail and work out a 
typical example. 



A Elaboration of the proof of the Theorem 4.2 



In this appendix, in Lemma A. 2 we prove (|4.5h ; in Lemma A.l we give the 



additional formulas beyond (4.16)-(4.21) necessary in proving the full Corollary 



Lemma A.l The Hirota symbols corresponding to the coefficients of Lemma 4-1 
with the noncontributing odd terms removed are, up to a constant, as follows 

Y 3 : -40^ + 30f + 0? 

Y 4 : -30104 + 20 2 3 + d 2 d\ 

Y 5 : |0 2 04 - |0i0 5 + |0I + ^0i 3 03 + |0i 2 0l + wA (A.l) 

Y 1>4 : -|0 2 4 + ]^0105 + ig0^ - 3^0? 03 - 3^01 

4Y lj4 + 10Y 5 : \0 2 d 4 - 20!0 5 + |0 3 2 + 10303 + \d\d\ 



whose action on r or yields the following differential equations for U = logr. 
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Y 3 
Y 4 

Y 5 



dp + Q(dfU) 2 + 3d%U - AdtdsU = 

-3<9i<9 4 £7 + 2d 2 d 3 U + dfd 2 U + 6(<9 2 l0(<9i<9 2 [/) = 

-^d&U + ±-%U + 6(dfUf + 3(dtU)(dfU) 
5 10 



+ 9d 2 d A U + AdlU + Adfd 3 U + 24(dfU)(d 1 d 3 U) 
+ 9{dfU)(d%U) + \d\d\XJ + 18(di<9 2 £/) 2 = 

Y M : —dxdsU + \d\U + 12(d 2 f/) 3 + 6(^7) (^£7) + 9<9 2 <9 4 f7 

5 5 

- 4<9f[7 + 2d\d 3 U + 12(dfU)(d 1 d 3 U) = 

4y m + ioy 5 : -Ad^u + d 2 d±u + plu + ld\d 3 u + Aidlu^d^u) 



+ S^fC/ + 4(01 d 2 ?7) 2 + 2(dfU)(diU) = 0. 



(A.2) 



Proof: Equations A.l follow immediately from the definition of Schur polynomials 



and (4.2), while (A.2) follows from (A.l ) and the definition of the Hirota symbol, i 



As a typical example of explicitly carrying out of the program given in the 



proof of Theorem 4.2 we prove (4.5) of Corollary 4.3 



Lemma A.2 For all p > 2, the Y 3 -KP hierarchy member yields the following 
equation for Q; 



<9 4 Q + 6(<9 2 Q) 2 + 5 2 , p (2-4£)dQ 



+(1 - 5 2 , P )(3«9 2 Q - (12(2=i)Vi0 + 4(1 - 5; 



Proof of Lemma A.ty From (A.2) and (4.16)-(4.21) conclude that g = lnre(t) 
and lnr(t) (remember r(t) = Tq(£) of Lemma 3.7) satisfy 
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Y 3 : dfg + Q(dfgY + 33&- 43^ = 0, 



d\g = d 2 g-(l-5 2 , p ) 



p — 1 
P 



tp-i, d{g = d 4 g, 



d.dsg = -(2e-l)dg, p = 2, 



did 3 g 



(dzdg - - 3)V_ 3 (1 - (1 - 6 3>p ), p > 2 



^ = d 2 2 g(l-5 2>p ), 
and substituting the last 4 relations into the first yields 



r 3 (<?) := Pg+G^Pg- ^ Vi(l-M) + 3(1 - S 2p )d 2 2 g 



-4(1 - 5 3p ) 

Set g = g = hir and then, since 9g = eg = 0, conclude that 



r 3 (<? ) = 6(i-(J 2|P ) 



p — l 



* P -iJ + 3(1 - h, P )d 2 go 
+12(1 - <5 3 , p )(l - <J 2>P ) ) t p _ 3 = 0. 



Since g = \nr E {t) = Q + lnr^^) := Q + and <9g = ego = 0, conclude from 
the above that 
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o = r 3 (Q + </o)-r 3 (</o) 



<9 4 Q + 6(9 



12(1 - <J. 



2,pj— ^p-i<9 2 (' 



+ 6 



p 
p 



Vl (1 - ^2,r 



+ 3(l-5 2iP )9 2 2 Q + 3(l-5 2iP )^ 



- 4(1 - 5 3 , p 



M e - + h, P )d 3 



+ 12(1 - 5 3jP )(l - fc,)^— ^V-s - T 3 (g ) 

p 

d 4 Q + Q(d 2 Q) 2 + ( -12^^Vi<9 2 Q + 3d 2 2 Q^ (1 - S 2 , p ) 



4(1 - ^ 



a? + (i - 5 2p )a 



which proves Lemma |A.2 



The rest of corollary 4.3 is proven in the same fashion using Lemma A.l and 



(4.16)-(4.21) 



References 

[1] M. Adler, M. Cafasso, and P. van Moerbeke. From the Pearcey to the Airy 
process. Electron. J. Probab., 16(36):1048-1064, 2011. 

[2] M. Adler, M. Cafasso, and P. van Moerbeke. Nonlinear PDEs for Fredholm 
determinants arising from string equations. arXiv:1201 '.6341- 

[3] M. Adler, T. Shiota, and P. van Moerbeke. A Lax representation for the 
vertex operator and the central extension. Comm. Math. Phys., 171(3):547- 
588, 1995. 

[4] M. Adler, T. Shiota, and P. van Moerbeke. Random matrices, Virasoro 
algebras, and noncommutative KP. Duke Math. J., 94(2):379-431, 1998. 



50 



[5] Mark Adler, Jonathan Delepine, and Pierre van Moerbeke. Dyson's nonin- 
tersecting Brownian motions with a few outliers. Comm. Pure Appl. Math., 
62(3):334-395, 2009. 

[6] Mark Adler, Jonathan Delepine, Pierre van Moerbeke, and Pol Vanhaecke. 
A PDE for non-intersecting Brownian motions and applications. Adv. Math., 
226(2):1715-1755, 2011. 

[7] Mark Adler, Nicolas Orantin, and Pierre van Moerbeke. Universality of the 
Pearcey process. Physica D, 239: (924-941), 2010. 

[8] Mark Adler and Pierre van Moerbeke. PDEs for the joint distributions of 
the Dyson, Airy and sine processes. Ann. Probab., 33(4):1326-1361, 2005. 

[9] Mark Adler and Pierre van Moerbeke. PDEs for the Gaussian ensemble with 
external source and the Pearcey distribution. Comm. Pure Appl. Math., 
60(9):1261-1292, 2007. 

[10] M. Bergere and B. Eynard. Universal scaling limits of matrix models, and 
(p,q) Liouville gravity. arXiv.0909.085J,.. 

[11] M. Bertola and M. Cafasso. Riemann-Hilbert approach to multi- 
time processes; the Airy and the Pearcey cases. Physica D, 
doi:10.1016/j.physd.2012.01.003, 2012. 

[12] M. Bertola and M. Cafasso. The Riemann-Hilbert approach to the tran- 
sition between the gap probabilities from the Pearcey to the Airy process. 
International Mathematics Research Notices, 2012(7):1519-1568, 2012. 

[13] P. M. Bleher and A. B. J. Kuijlaars. Random matrices with external source 
and multiple orthogonal polynomials. Int. Math. Res. Not., (3): 109-129, 
2004. 

[14] A. Borodin and P. Deift. Fredholm determinants, Jimbo-Miwa-Ueno r- 
functions, and representation theory. Comm. Pure Appl. Math., 55 (9): 1160- 
1230, 2002. 

[15] E. Brezin and S. Hikami. Correlations of nearby levels induced by a random 
potential. Nuclear Phys. B, 479(3) :697-706, 1996. 

[16] E. Brezin and S. Hikami. Extension of level spacing universality. Phys Rev., 
E 56:264-269, 1997. 

[17] E. Brezin and S. Hikami. Level spacing of random matrices in an external 
source. Phys. Rev. E (3), 58(6, part A):7176-7185, 1998. 



51 



[18] E. Brezin and S. Hikami. Universal singularity at the closure of a gap in a 
random matrix theory. Phys. Rev. E (3), 57(4):4140-4149, 1998. 

[19] T. Claeys, A. Its, and I. Krasovsky. Higher-order analogues of the Tracy- 
Widom distribution and the Painleve II hierarchy. Comm. Pure Appl. Math., 
63(3):362-412, 2010. 

[20] Etsuro Date, Masaki Kashiwara, Michio Jimbo, and Tetsuji Miwa. Trans- 
formation groups for soliton equations. In Nonlinear integrable systems- 
classical theory and quantum theory (Kyoto, 1981), pages 39-119. World 
Sci. Publishing, Singapore, 1983. 

[21] Etsuro Date, Masaki Kashiwara, and Tetsuji Miwa. Transformation groups 
for soliton equations. II. Vertex operators and r functions. Proc. Japan Acad. 
Ser. A Math. Sci, 57(8):387-392, 1981. 

[22] J. Harnad and A. R. Its. Integrable Fredholm operators and dual isomon- 
odromic deformations. Comm. Math. Phys., 226(3) :497-530, 2002. 

[23] G. F. Helminck and J. W. van de Leur. An analytic description of the vector 
constrained KP hierarchy. Comm. Math. Phys., 193(3) :627-641, 1998. 

[24] A. R. Its, A. G. Izergin, V. E. Korepin, and N. A. Slavnov. Differential 
equations for quantum correlation functions. In Proceedings of the Confer- 
ence on Yang-Baxter Equations, Conformal Invariance and Integrability in 
Statistical Mechanics and Field Theory, volume 4, pages 1003-1037, 1990. 

[25] Kurt Johansson. Discrete polynuclear growth and determinantal processes. 
Comm. Math. Phys., 242(l-2):277-329, 2003. 

[26] Kurt Johansson. The arctic circle boundary and the Airy process. Ann. 
Probab., 33(l):l-30, 2005. 

[27] V. Kac and A. Schwarz. Geometric interpretation of the partition function 
of 2D gravity. Phys. Lett. B, 257(3-4) :329-334, 1991. 

[28] Igor Moiseevich Krichever. General rational reductions of the Kadomtsev- 
Petviashvili hierarchy and their symmetries. Funktsional. Anal, i Prilozhen., 
29(2):l-8, 95, 1995. 

[29] T. Miwa, M. Jimbo, and E. Date. Solitons, volume 135 of Cambridge Tracts 
in Mathematics. Cambridge University Press, Cambridge, 2000. Differential 
equations, symmetries and infinite-dimensional algebras, Translated from the 
1993 Japanese original by Miles Reid. 



52 



[30] A. Okounkov and N. Reshetikhin. On the birth of a random process. Mosc. 
Math. J., 6:553-566, 2006. 

[31] Michael Prahofer and Herbert Spohn. Scale invariance of the PNG droplet 
and the Airy process. J. Statist. Phys., 108(5-6): 1071-1106, 2002. Dedicated 
to David Ruelle and Yasha Sinai on the occasion of their 65th birthdays. 

[32] Igor Rumanov. All the lowest order pde for spectral gaps of gaussian matri- 
ces. arXiv: 1008.3560. 

[33] Igor Rumanov. Universal Structure and Universal PDE for Unitary Ensem- 
bles. J. Math. Phys., 8, 8 2010. 

[34] Graeme Segal and George Wilson. Loop groups and equations of KdV type. 
Inst. Hautes Etudes Sci. Publ. Math., (61):5-65, 1985. 

[35] Craig A. Tracy and Harold Widom. A system of differential equations for 
the Airy process. Electron. Comm. Probab., 8:93-98 (electronic), 2003. 

[36] Craig A. Tracy and Harold Widom. Differential equations for Dyson pro- 
cesses. Comm. Math. Phys., 252(1-3) :7-41, 2004. 

[37] Craig A. Tracy and Harold Widom. The Pearcey process. Comm. Math. 
Phys., 263(2):381-400, 2006. 

[38] Johan van de Leur. The vector A;-constrained KP hierarchy and Sato's Grass- 
mannian. J. Geom. Phys., 23(l):83-96, 1997. 

[39] Wolfgang Wasow. Asymptotic expansions for ordinary differential equations. 
Dover Publications Inc., New York, 1987. Reprint of the 1976 edition. 

[40] Harold Widom. On asymptotics for the Airy process. J. Statist. Phys., 
115(3-4):1129-1134, 2004. 



53 



